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ABSTRACT 

A migrating planet can capture planetesimals into mean motion resonances. However, resonant 
trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple 
Hamiltonian system for first and second order resonances, we explore how the capture probability 
depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling 
factors as a function of the planet mass and resonance strength to estimate the planetary migration 
rate above which the capture probability drops to less than 1/2. Applying our framework to multiple 
extra solar planetary systems that have two planets locked in resonance, we estimate lower limits for 
the outer planet's migration rate allowing resonance capture of the inner planet. 

Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation 
subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We 
present factors that can be used to estimate this critical planet eccentricity. Applying our framework 
to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would 
make its 2:1 resonance fail to capture twotinos. The capture probability is affected by the separation 
between resonant subterms and so is also a function of the precession rates of the longitudes of periapse 
of both planet and particle near resonance. 
Subject headings: celestial mechanics 



1. INTRODUCTION 

Resonances can capture particles in slowly varying 
dynamical systems. For example, a planet migrating out- 
ward can trap planetesimals in resonances exterior to it; 
as Neptune trapped the Plutinos in the Kuiper bel t (e.g 
P^rnandez & Id 1984; Malhotra 1995; Hahn & Malhotra 
1999tilda et al.>2000ciCh iang & Jordan 2002; Zhou et al 



20021 iLevison fc Morbid elh 2003; Wyatt 2003:). A planet 



migrating inward can trap planetesimals or planets in 
resonances interi or to it (e.g., Ouillen & Holman 2000; 
iKlev et al] l2004|) . Dust spiraling inward under dissi- 
pational forces can become trapped in exterior mean 
motio n resonances with a planet (e.g.. iSicardv et aT 
19931 iDermott et all [l99l Marzari fc Vanzini 
1994 ; iGou fc Zook 119971 11999 : Ozernov et jd 

t 



20001 \m\u^eA'^?2m' 
Deller fc Maddisonll2005|)T 



iMoro-Martin 
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An elegant and predictive theory of resonant capture 
has been developed for adiabatically varying non-chaotic 
integrable resonant systern s (Yoder 1979; Henrard 1982; 
iHenrard fc LemaitrdllEisl lMalhotralll99(]|) . and for the 
same systems b ut varying with weak nonconservative 
forces l)Gomesll 997'). This theory was first applied to sys- 
tems of tidally locked satellites fBorderies fc GoldreicW 
[1984; Pcalc i986; Dermott et al. 1988; Malhotra 1990). 
However this theory does not apply to systems that 
are near or in the non-adiabatic regime, or are chaotic. 
Numerical explorations of drifting and migrating sys- 
tems have revealed differences between measured capture 
probabilities and those predicted by the adiabatic theory. 
For dust drifting under dissipational forces, small parti- 
cles can be drifting sufficiently fast that they are in the 
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non-adiabatic regime. In thi s case, the capture proba- 
bility is reduced (e.g., [Gomes 1995; .Liou fc Zook 1999). 
Numerical simulations of Neptune's migration show that 
if Neptune migrates ra pidly, the capture prob ability of 
resonances is red uced |ld"a et al.l[2000t iFriedla nd 200 3 
IChia,ng fc ■T^^7d^l200l . 

In dynamical systems chaotic motion can arise from 
the overlap of resonances (e.g., Wisdom 1980:!). Mean 
motion resonances associated with perturbations by a 
planet are ma de up of more than one resonant pertur- 
bative term. Holm an fc Murravl l|1996j) showed how to 
predict the Lyapunov times of asteroids from the over- 
lap of these resonant terms. The complex behavior of 
chaotic resonances could influence the way they cap- 
ture particles when they are varyin g (e.g. .IDermott et alJ 
r^88; 'Tittemore & Wisdom Iggo"; "Sicar dv et all 119931 
[^r zari fc Vanzini 1994; Ouillen 2001). For exam- 
ple, iTfttemore fci'Wisdoml ljl99Ql found that the drift- 
ing chaotic resonances of the Uranian satellites exhibited 
different behavior than non-chaotic resonances. 

In this paper we strive to develop a general frame- 
work that will allow us to better predict the capture 
probability of resonances. We would like to know when 
they are likely to capture and how the capture prob- 
ability depends on the drift or migration rate, reso- 
nance order and strength, initial particle orbital ele- 
ments, and size and separation of subresonances. Some 
of these parameters depend on the planetary proper- 
ties and others depend on the particle properties. We 
would like a flexible framework that could be used to 
place constraints on migrating planetary extrasolar sys- 
tems and the properties of the planets and planetesi- 
mals residing within them. Previous works have shown 
how to predict the capture probability as a function of 
initial particle eccentricity in the adiabatic limit (e.g.. 



I Henrardl ITosl IBorderies fc CToldreichl ITosl IMalhotral 
Il99nt iMurrav & Dermottlll999D . Here we strive to un- 
derstand two regimes that have not been well explored 
for mean motion resonances. We probe the nonadiabatic 
regime in which the drift rate is so fast that the reso- 
nance can fail to capture. We also explore resonances 
comprised of multiple resonant subterms. 

Our approach is to understand the simplest Hamil- 
tonian model that can illustrate the dynamics of reso- 
nance capture. In section 2 we formulate the problem 
in terms of an idealized Hamiltonian with multiple res- 
onant terms. By rescaling, this simple model allows dif- 
ferent resonant systems to be treated in a similar way. 
In section 3 we consider the system with one resonant 
term and measure the capture probability as a function 
of drift rate and initial particle momentum. We extend 
previous analytical work bylFriedland (20^ for the first 
order resonances and find general expressions that ac- 
count fo£_trends exhibited by previous simulations (e.g., 
IWvatti I2003D . In section 4 we measure capture proba- 
bilities for first and second order resonances containing 
multiple resonant terms. In section 5 we illustrate how 
our framework can be applied to planetary systems. A 
summary and conclusion follows. The appendix lists co- 
efficients for common resonances so that drift rates can 
be predicted and compared to numerical studies. 

2. HAMILTONIAN FORMULATION 

We employ the Poincare coordinates 

X = M + w, 7 = —Tu 

and their associated momenta 

L = ^GM^a, r = y/GAha{l - ^1 - e^) 

where Af, is the mass of the star, A is the mean longitude, 
M is the mean anomaly, w is the longitude of pericenter, 
a is the semimajor axis, and e is the eccentricity. These 
variables are those describing the orbit of a particle or 
planetesimal in a plane. The Hamiltonian for the Kep- 
lerian system in these coordinates restricted to a plane 
is 

H{L,\,V,^) = -^-R 

where R is the disturbing function that depends on the 
coordinates of the particle and on the coordinates of the 
planet. The planet's semimajor axis and mass are de- 
noted ap and m^, respectively. The planet's other co- 
ordinates are subscripted in the same way. The mean 
motion of the particle n — \ where A is the derivative 
with respect to time of A. 

Hereafter we adopt a unit convention with distances 
in units of the planet's semi-major axis, ap, at a time 



to. Time is put in units of y^a^/GM*. We define /i 

to be the mass ratio = rUp/M^. At low eccentricity, 
T/L « e^/2, relating the momentum F to the particle 
eccentricity. We often give the particle semimajor axis 
in terms of the variable a = ^ \i a > ap (external to the 
planet) and a = ^ for a < (internal to the planet). 

The unperturbed Hamiltonian or that lacking the dis- 
turbing function 

i7o(i,A;r,7) = -^ 



We consider the the j '■ j — k exterior mean motion res- 
onance (planet is an interior perturber). We perform a 
canonical transformation using the mixed variable gen- 
erating function 

F2 - I{jX - U - k)Xp) 
leading to new variables 



and new Hamiltonian 



V' = jA - (j - fc)Ap 



i7^(/,V;r,7) = 



-1 

2pP 



- U - k)h 



,p. 



We now expand around the resonance. Let 

A = / - /o 

and 

1 

{j - k)np{to). 



(1) 



■J -'o 

Since we have adopted units np{tQ) = 1, we find Iq — 

01-^^ I j where a = ^ = ('^-j-j on resonance, as 
expected. Our Hamiltonian now reads 



Xo(A, t/i; r, 7) = constant — [j — k){np — 1)A - 
We can write the unperturbed Hamiltonian as 



3A2 



ifo(A,7^;r,7) 
with coefficients 



oA -|- 6A -I- constant 



3-2 2 

a~ 7 a 

b^-ij~k){np-l). (2) 

We now recover the disturbing function that is is tra- 
ditionally expanded as a cosine series of angles in or- 
ders of planet and particle eccentricity. We keep the 
terms inducing precession of the longitude of periapse 
and low order terms (in eccentricity) containing ip. The 
full Hamiltonian 



K{A, iP; r, 7) = aA^ + 6A + cF 



(3) 



■ J2 '^fe.pr^'^"^^/' cos (^ - (fc - p)w - pwp) 



p=0 



with coefficient 



-/i2/2a 



1/2 



(4) 

2V/L ~ 2Ta^/^. 



We have used the approximation 
Here the perturbation strengt hs, 5k..n, are functions of a , 
j and Laplace coefficients fsee lMurrav fc D ermott'1999'). 
The c term describes secular precession of the longitude 
of periapse and depends on the functio n /2 given in the 
appendix bv IMurrav fc DermottI ljl999|) and is evaluated 
at a with index j = 0. As have previous studies, we 
have neglected the dependence of a on time as the planet 
migrates, and cosine terms from the disturbing function 
expansion whic h are expected to average to zero near 
resonance (e.g. .IBor deries fc Goldreichlll984tlPealdll986t 
Holman fc Murray 199(i HMurrav fc Dermottlll999D . The 
perturbation strengths depend on the planet mass and 
eccentricity as 5k,p <x. pe^. More detailed expressions are 
listed in the appendix. 



The a bove canonical tran sformations are similar to 
those of iHolman fc Murravl ifT? 



except we have fo- 
cused on resonances exterior to a planet rather than those 
interior to it. We have also explicitly kept the 6A term. 
A system with a migrating planet would be described 
by a time dependent b coefficient. This allows us to ex- 
plore the dynamical behavior as particles pass through 
resonances . 

Holman fc Murravl l|1996f) showed that the above 
Hamiltonian is similar to a periodically forced pendu- 
lum and that the overlap of the different resonant terms 
of Equation Q can induce largescale chaotic behavior. 
However most previous explorations of resonant capture 
have only considered one dominant resonant term. When 
the migration rate is slow, the adiabatic theory developed 
by (Yodcr 1979 .: Hcnrard 1982; Bordcrics & GoldrciclJ 
[K)84; Mal hotralTl99 0: Murrav & Dcrmott 19990 applies. 
In the next section we explore this simpler situation, but 
allow the migration rate to be fast or non-adiabatic. 

3. PROBABILITY OF CAPTURE IN A SINGLE 
RESONANCE 

In this section we explore the simpler Hamiltonian con- 
taining only one dominant resonant term. This situation 
would be appropriate if the planet's eccentricity is very 
small in which case the 6kfi term dominates. The Hamil- 
tonian (Equation 121) including only this term 

K{A,TP;T,j)^aA^ + bA + cT (5) 

+4,or''/^ cos {tp - kw) 

It is convenient to perform a canonical transformation 
with generating function 



leading to new variables 

Ji = r, 

and new Hamiltonian 
K'{T, 



t 
' k 



-^2» = a ( ^ 



£2 

If 

2aJ2 



r + 6J2 



+4,or'-'/' cos(fc0) 



(6) 



Note that J2 is conserved and is small for initial con- 
ditions near resonance with small initial eccentricity (or 

r). 

Dropping constant terms and setting b' — (2a J2 + b) + 
kc, the Hamiltonian in Equation ((BJ 



+4,or^/'cos(fc0). 

By rescaling momentum and time 



-2/(4-fc) 



^ 14,0 



2/(4-fe) 



r 

(2-fc)/(4-fc) 



(7) 



we can write this as 

Kit, = f 2 + &f + (-l)'=f '^/^ cos(fc0) (8) 

where 

(fe-2)/(4-;e) 



6 = 6'|4,or2/(4-'=) 



sign(a). (9) 



We relate the drift rate b from the migrating planet to 
that of our scale-free system 



db 




b 






k 



|4,o| 



-4/(4-A;) 



2(fe-2)/(4-fe) 



(10) 



The form of the Hamiltonian in Equation Q is 
(excepting for factors of \/2)^ identical to that 
used to explo r e capture in the adiabatic li mit (e.g., 
HenrardI 119821: iBorderies fc GoldreichI 119841: iMalhotral 
1990HMurrav fc Dermottlll999j) . 



3.1. Capture probability as a function of drift rate and 
initial momentum for first and second order 
resonances 

We ask the following question: Above what drift rate 
(^) does the resonance fail to capture ? The system be- 
haves adiabatically when it takes longer than an oscilla- 
tion period for the system to pass through the resonance. 
For low initial momentum, the width of the resonance is 
~ 1, and the period of oscillation is ~ 1. Consequently 
we expect the system would evolve adiabatically when 

^1. To go beyond this limit and estimate the 

probability of capture as a function of drift rate in the 
non-adiabatic regime, we would need to find solutions to 
Hamilton's equation. It is non-trivial to find solutions 
to the equations of motion for Equation ((SJ when b is 
a function of t ime and the sys tem is not varying adia- 
batically (e.g., iFriedlandl 120011) . Consequently we have 
integrated Hamilton's equations of motion numerically 
to explore the non-adiabatic regime. Once we numeri- 
cally understand the behavior of the scale-free Hamilto- 
nian (Equation |HJ|, we can make predictions for systems 
in the same form using the factors of Equation ^ . 

Our procedure for numerical integration is as follows. 
Hamilton's equations for Equation JHJ are integrated us- 
ing a conventional Burlisch-Stoer numerical scheme. The 
initial angle is chosen randomly. We assume that b is pro- 
portional to time so that only one parameter ^ specifies 
the time dependence of the system. To ensure that the 
particles were initially outside of resonance we require 
the initial value of |6| to exceed 1. The parameter b was 
initially chosen to be ~ —15, well outside the resonance. 
The system passes through resonance when 5 ~ so the 
timescale until capture is t capture ~ &mit( J^)"^- The sys- 
tem is integrated at least twice the capture time. Two 
sample integrations are shown in Figure ^ Figure ^ 
shows an integration illustrating a particle that is cap- 
tured into resonance and Figure^ shows one where no 
capture takes place. 

In the adiabatic limit, the capture probability is 1 when 
the initial particle eccentricity is smaller than a limiting 

^ By rescaling our momentum by a factor of 2*^+^/'''"*') , the 
Hamiltonian becomes + bV + (-i)k2(k+2)/2f^k/2 cos{k(f>) as by 
IMurrav fc DermottI 119991) . 



4 



value, eiim, that depends on the order and width of the 
resonance Ce.g iHmraTdl llQSl IBorderies fc ColdreichI 
IT981 iMalhotTa 1199(11 IMurrav & DermottI 119991) . How- 
ever when the drift is not adiabatic, the capture prob- 
abihty could depend on the initial particle eccentricity 
(or momentum F) even when it is below this limiting 
value. Consequently we measured capture probability 
as a function of both drift rate (^) and initial parti- 
cle momentum. For each value of drift rate and initial 
momentum, we integrated the system 100 times (each 
time with a different randomly chosen angle) to measure 
a capture probability. 

After the system captures, the momentum increases 
with time (see Figure ^) and the resonant angle 4> li- 
brates about a fixed value (0 or tt) rather than circulat- 
ing. If no capture takes place, the momentum jumps as 
the system pass through resonance (see Fig. When 
the momentum increases to a value exceeding the reso- 
nance width (approximately 1 as we have rescaled the 
Hamiltonian) we identify the system as having captured 
into resonance. We used a limiting value of F = 5 to 
identify captures. In resonance the angle librates around 
a fixed value. The condition cj) ^ implies that ^ ~ 0, 
and f ~ —6/2 in resonance. We use this condition 
to ensure that we integrate Hamilton's equations long 
enough that the momentum crosses our limiting momen- 
tum value when the resonance captures. 

Figure 13 and 13 show capture probabilities that we have 
measured numerically for first and second order reso- 
nances {k = 1 and 2). In the adiabatic limit for T{to) 
below 3/2 for fc = 1, or 1/8 for fc = 2, the capture prob- 
ability is one.^ Above these limiting initial momentum 
values the capture probability is less than 1 in the adia- 
batic limit. For f (to) = 2.3,_shown as stars in Figure[21 
above the limiting value of Fo^Hm = 1-5, we see that the 
capture probability never reaches 1. At lower drift rates 
the capture probability approaches a constant value for 
this initial condition, consistent with the prediction in 
the adiabatic limit. The same behavior is seen for fc = 2 
with initial momentum T{to) = 1 (shown as stars in Fig- 
ure ISJ. This momentum is eight times the limiting value 
which ensures capture in the adiabatic limit. 

To quantify the width of the probability function, we 
have fit a function to the capture probability p{u) — 

0.5 (^1 - tanh (^ "~^^^^ ^^ as a function of u = log^Q |^|. 
Here is the log of the drift rate at which the cap- 
ture probability is 1/2, and w describes the width of the 
drop. For large w, the slope is shallow, for small w the 
drop is a steep function of the drift rate. The drift rates 
at which the capture probability is half and a quarter, 
and the widths of the probability functions are shown 
as a function of initial momentum in Figures ^ and |5l 
for first and second order resonances, respectively. For 
initial momentum sufficiently low, (e.g., 10~^ for fc = 1 
and 10^^ for k — 2) the drift rate at which the proba- 
bility is half approaches a limiting value. The steepness 



^ O ur momentum is half or 1/4 times that of 
119991 ) for fc = 1 and fc = 2 respectively. 



Murrav K 


z, Dermott 


Murray & 


L Dermoti 



I199HI) list critical momentum values of 3 and 1/2. These crit- 
ical values are used to find the maximum parti cle eccentricity, 
eHm,, that ensures capture i n the adiabatic limit IHenrar dl 11981 
IBorderies fc GoldreichllT98l : IMalhotrallT99nh . 



of the transition between 100% capture and 0% capture 
is narrower in its range of drift rates at lower values of 
initial momentum. The lower the initial momentum, the 
sharper the transition between a capture probability of 1 
and zero. A sharp transition is reached at a lower initial 
momentum for fc = 2 than in the fc = 1 case. For initial 
momentum near Fo^zim, the limiting value ensuring cap- 
ture in the adiabatic limit, there is a regime or a range of 
drift rates where the_ capture probability is intermediate. 
In other words, for r(io) ^ 1 the widths w ~ 1. 

We now consider the situation where the transition be- 
tween a probabilitj^ of 1 and is sharp. This is true for 
initial momentum F(to) ^ 10~^ and 10~^ for fc = 1 and 
for fc = 2 respectively. For these initial momenta we 
measure the critical drift rate where the transition takes 
place. From our numerical integrations, the dynamical 
system fails to capture for drift rates faster than 



2.0 



'0.25 



for fc = 1 
for fc = 2. 



(11) 



For first order resonances the drift rate for a capture 
probability of 1/2 is not strongly dependent on the ini- 
tial momentum as long as this lies below Fo.iim- However 
for second order resonances when the initial momentum 
r(to) ~ 1, the capture probability of half occurs at a 
drift rate that is about 10 times that at low initial mo- 
mentum. We have approximated the dependence of the 
half probability drift rate on the initial momentum with 
the following function (shown as a dotted line in Figure 

13 



d6 

dT 



0.25 1 



1/2 



f(io) 



3 X 10- 



0.25 



(12) 



The power, 0.25, is not necessarily theoretically meaning- 
ful. This function is a reasonable match to the measured 
points of Figure |S1 for initial momentum r(to) < 1. For 
higher initial momentum, r(to) ^ 4, the probability of 
capture never exceeds 1/2. This initial momentum (4) 
exceeds f o,;™ by a factor of 32. 

In the adiabatic limit, the probability of capture drops 
as a function of increasing momentum (or eccentricity) 
when the initial momentum is above VQMm- However the 
probability drops faster for fir st order re sonances than 
for second order resonances (iHahn fc Malh otra 1999). In 
the adiabatic limit, the probability of capture for a first 
order resonance drops to 1/2 for f (tg) ~ 2.3 (less than 
twice the limiting value of 1.5), whereas for second order 
resonances the probability of capture drops to 1/2 for 
F(to) ^ 4 or 32 times the limiting value. 

Using the critical drift rates for ^ (equation lll|l we 
can invert Equation (|10|l to determine which resonances 
can capture at a particular drift rate. We find that reso- 
nances are likely to capture for h slower than the critical 
rates 

~2|<5i^o|'/=^|ap/' forfc = l 

0.5(5|_o for fc = 2. (13) 

Using Equation ||2Jl to replace h with the planet's mean 
motion 

|np,c„t|^2(j-l)|5i,o|4/3|„|2/3 forfc = l 

0.5(j - 2)(5|o for fc = 2. (14) 



5 



For second order resonances, the rate given above can 
be modified by the function given in Equation l|12|) to 
estimate the rate at which the probabihty is 1/2 as a 
function of initial momentum; 

iwi-o.5o--2)4,(i+ j;;;t^^j °-" (15) 

where eq is the initial particle eccentricity. This is ex- 
pression is valid for initial particle eccentricities smaller 
than Co ^ lOejim where eum is the eccentricity limit en- 
suring capture in the adiabatic limit. An expression for 
eiim is given in the appendix. The factor of 10 comes 
from the range covered by the function shown in Equa- 
tion (|12|1 . The curve shown in Figure [S] is a reasonable 
match up to r(to) ~ 1 which is approximately lOTo^um- 

The above relations fEauationsll4l and ll5|) allow us to 
estimate the likelihood of resonance capture in different 
astronomical settings. The strength of the perturbative 
cosine terms iSk,p) are proportional to the planet's mass 
or /i, however the critical drift speed depends on /i^/"^ for 
k = 1 and on /x^ for /c = 2. We see that the critical drift 
rates for capture are strong functions of the planet mass 
and this is particularly true for the second order reso- 
nances. Slower drift rates are required to allow resonant 
capture for lower mass planets. 

The dependence of critical drift rate on planet mass 
provides a qualitative explanation for some features of 
numerical simulations which start with particles in ini- 
tially low eccentricity orbits. The above relation pre- 
dicts that only more slowly drifting particles will be able 
to capture into higher order resonances. We can under- 
stand why the 5:3 resonance requires slower migration 
rates than the 2:1 and 3:2 re sonances to capture in the 
simulations of iWvattI |2003). Although we have inte- 
grated a time dependent Hamiltonian system, we can 
expect similarities between this system and the slowly 
drifting non-conservative systems. Simulations of dust 
drifting inward via dissipative forces tend to show that 
large dust particles are captured into higher order res- 
onances than smaller pa rticles (e.g., Marzari & Vanzini 
ITflflllLion fc ZookllToM) . This follows since small par- 
ticles drift faster than larger ones and the higher order 
resonances require slower drift rates to capture. 

Based on results of simulations of Neptune's migration 
^a et al. (2000) proposed that the critical drift rate de- 
pended on planet mass to the 4/3 power. They restricted 
their study to fc = 1 resonances so their prediction is con- 
sistent with our previous equation. This po wer depen- 
dence was confirmed with analytical work by iFriedlandl 
II2001D. als o for the fc = 1 resonance, and numerical work 
bv IWvatH l|2003() . We confirm the steeper dependence 
on planet mass of the 5:3 resonance cap ture probability 
measured numerically bv iWvat^ l)2003j) and specifically 
predict that the critical drift rate is cx /x^ for second or- 
der resonances. The relation for the critical drift rate 
(Equation is both consistent with and more general 
than the scaling found by these previous studies. Be- 
cause we have related the critical drift rate (via scaling) 
to the resonance strengths, the formulation given here 
can be applied to any first or second order mean motion 
resonance. 

Here we have also found that the probability of capture 
when the drift is not adiabatic is a non-trivial function 



of initial particle eccentricity. The transition between 
a probability of 1 and becomes smoother (covering 
a larger range of drift rates) as the initial momentum 
approaches the minimum value ensuring capture in the 
adiabatic limit (see Figures El El EI and EJ • For first or- 
der resonances, the midpoint drift rate (corresponding 
to a probabihty of capture of 1/2) does not significantly 
depend on the initial particle eccentricity. However for 
second order resonances the midpoint is at a higher drift 
rate when the initial momentum or particle eccentricity 
is higher. The increase in drift rate with initial particle 
eccentricity allowing capture for second order res o nance s 
was described previously by iHahn fc Malhotral l)2005D . 
However we do not predict the same dependence on res- 
onance width and planet mass. This is because we have 
restricted our study to initial particle eccentricity near 
or below eum , and iHahn fc Malhotraj l)2005|) considered 
initial particle eccentricity exceeding eum- 

Our estimate for the critical drift rate above (Equation 
114(1 is appropriate for a wide range of initial particle ec- 
centricities for first order resonances (as long as they are 
below the limiting value, eum)- The half probability drift 
rate's dependence on the initial particle eccentricity can 
be estimated for second order resonances using Equation 
H15I) when the initial particle eccentricity is lower than 
~ 10 times eum- The framework we provide here can we 
used to estimate the half probability drift rate for any 
second order resonance. 

Numerical studies report intermediate probabili- 
ties for ca pture i nto first order resonances from 
simulations lllda et a l. 2000; Quillcn & Holman 200^ 
iChiang fc Jordanll2002l: IWvattll2003() . By intermediate, 
we mean not close to zero or 1, or at a ~ 50% level. 
Here we have found that the dynamical system described 
by Equation (jHJ for fc = 1 with only a single resonance 
term has a limited range of drift rates where the capture 
probability^ is intermediate, unless the initial particle mo- 
mentum r is of order 1. This regime corresponds to an 
initial particle eccentricity within a factor of a few of 
eiim, the limiting value ensuring capture in the adiabatic 
limit. The limiting eccentricity depends on the resonance 
strength, and planet mass to the power fc/(4 — fc) (using 
the squa re root of the s cale-free momentum in Equa- 
tion O; (IMalhotralllQQOtlMurrav fc Dermotl ll999'). For 
weaker resonances, the initial particle eccentricity limit is 
more restrictive. It is possible that some of the numerical 
simulations are effectively in the regime of intermediate 
capture for certain resonances due to their initial parti- 
cle eccentricity distribution. We return to this issue in 
later sections as we identify other regimes of intermediate 
capture probability for first order resonances. 

Though the limiting eccentricity is smaller for second 
order resonances, the probability of capture drops more 
slowly in the adiabatic limit as a function of initial parti- 
cle eccentricity at values above eum- Furthermore since 
the half probability drift rate increases with initial parti- 
cle eccentricity (Equation 1151) . higher eccentricity parti- 
cles can be captured at higher drift r ates than lower ec- 
centr icity particles (as pointed out bv IHahn fc Malhotral 
12005(1 . This also implies that the second order resonances 
have a larger regime in both range of initial eccentricity 
and drift rate where the probability of capture is inter- 
mediate. 
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4. THE ROLE OF AN ADDITIONAL RESONANCE TERM 

We now consider the differences in tfie dynamics of 
capture caused by the addition of a secondary resonant 
term for a first order resonance. We rescale the momenta 
and time for Hamiltonian given in Equation ljAl|l for k = 
1 according to the Equations Q. This gives us unitless 
momenta and time 

Ki(A>;f,7) = A2 + 5A + cf (16) 

— f "^/^ cos('(/' — tn) + ecos(V' — Wp) 

where 



i<5i,iii^i.or'/'kr/' 



'-c\Si.o\ 



-2/3|„|-l/3 



sign(a). 



(17) 



The coefhcients that can be adjusted are the width of 
the second resonance compared to the first (set by e) 
and the separation between the two resonances (set by 
c). Variation in Wp can be absorbed into our coefficient 
c. The term proportional to cos('0 — cc) is often called 
the e-resonance since T^^"^ (x e. The other term can be 
called an e'-resonance or a corotation resonance since it 
does not depend on the particle's longitude of perihelion 
or vu. 

_ Because the corotation resonance does not depend on 
r, it does not grow in volume as the planet migrates. Had 
we allowed a to depend upon time, the resonance width 
would grow slightly but not significantly as the planet 
migrates. Because the resonance volume in phase space 
does not grow as the planet migrates, thi s resonance 
should not capture particles ijYodeJ Il979() . However 
when this resonance overlaps the ot her, the system can 
exhibit large scale chaotic behavior fHolman fc Murravl 
IT99&) . Hence the coupling of the two resonant terms may 
influence the probability of capture into the e-resonance. 
We now ask: what is the capture probability of the 

above Hamiltonian as a function of drift rate, and 
secondary perturbation strength, e ? To answer this ques- 
tion we numerically integrate Equation Ijltill , for different 

parameters ^ , and e. In this section we work in the limit 
of low initial particle momentum. This ensures that the 
probability of capture is 1 in the adiabatic limit and that 
the transition between capture and no capture would be 
a sharp transition of drift rate in the case of a single 
resonance. 

Our procedure for numerical integration is the same as 
described in the previous section. Initial angles are ran- 
domly chosen. Initial momenta (A, F) are set to small 
values to ensure a sharp transition when e = and the 
initial momentum are small. The parameter b is ini- 
tially chosen to be ~ —15. For each value of ^ and e 
we integrated the system 100 times to estimate a reso- 
nant capture probability. After the system captures, the 
momenta variables increase with time and the resonant 
angle (p = ijj — w librates about a fixed value. If no cap- 
ture takes place, the momenta jump as the system pass 
through resonance. 

The coupled two dimensional system exhibits different 
dynamics than the one-dimensional system considered 
in the previous section. For example, the resonance can 
capture for a short period of time, a trajectory we refer 
to as a temporary capture. An example of a simulation 
that illustrates a temporary capture is shown in Figure 



We find that temporary capt ures tend to oc cur for 
larger values of e and drift rate. iQuillenI l|20niD previ- 
ously showed that temporary capture was exhibited by 
overlapped resonances using a similar drifting Hamilto- 
nian model. 

During a temporary capture, the momenta increase. 
This also happens if the particle is captured. As in the 
previous section, we identify a capture if the momentum 
at the end of the integration exceeds a value of 5. How- 
ever we then reclassify the integration as a temporary 
capture if the momenta lie below that expected from a 
particle still in resonance. These two situations can be 
differentiated because a particle in resonance has mo- 
mentum proportional to the time since capture. Tem- 
porary captures are excluded when we calculate the cap- 
ture probability. However, had we integrated the systems 
longer, it is possible that a particle identified on a short 
timescale as captured would later drop out of resonance. 
In other words the precise fraction of captures is depen- 
dent on the timescale over which we have integrated these 
systems and the value of momentum that we have used 
as a limit to identify captures. This makes our capture 
probability numerical measurements uncertain primarily 
at high values of drift rate \ > 1 , and large secondary 
perturbation strength e ^ 1, the regime where we have 
found temporary captures to be more common. 

In Fig. Owe show a contour plot of the resonant cap- 
ture probability for k — 1, initial T{to) = 10"'' and res- 
onance separation c — 0. For high values of the drift 
rate and low values of the secondary perturbation, e, (on 
lower right in these contour plots) the transition between 
capture and no capture is sharp and happens at the crit- 
ical drift rate measured in the previous section. However 
for lower drift rates and higher values of e the resonance 
fails to capture for e > 1. For low drift rates, the tran- 
sition between capture and no capture is also sharp, but 
is a function of e instead of drift rate. There are two 
regimes, that where the drift is so fast that it fails to 
capture, and that where the corotation resonance is so 
large that it prevents capture. 

For first order resonances we consider the possibility 
that e could be of order 1. Since e cx |5i,o|~'*^'^, |^i,o| « 
fi and |(5i,o| c>c /iCp, the coefficient e (x /i~^/'^ep. For 
small planet mass /i we see that the rescaled secondary 
perturbation strength could be high even at moderate 
planet eccentricity Cp. The coefficient e could be of order 

1 particularly if the planet eccentricity is moderate. 
Why is it that for e > 1 the e-resonance fails to capture 

particles? A possible explanation is that an increase in F 
caused by the corotation reduces the capture probability, 
in the same way that an increase in the initial momen- 
tum value does. We would expect that an increase of 
initial momentum of size ~ 1 caused by the corotation 
resonance when e ~ 1 would strongly reduce the capture 
probability. This is consistent with the limiting value 
of e < 1 for capture. This qualitative explanation is 
also consistent with the lack of dependence of the criti- 
cal value of e on drift rate (see Figure O at drift rates 

^ < 1. We suspect that the corotation resonance pre- 
vents capture into the e-resonance because the corota- 
tion resonance raises the particle eccentricity during the 
resonance encounter. If this were true, then we would 
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expect that the capture probabiUty would be influenced 
by the resonance separation. Up to this point we have 
only considered resonances with separation c = 0. 

4.0.1. Separated first order resonances 

To further explore the role of multiple resonant terms, 
we consider the situation when the resonance separation 
is non zero (c not small). The order that the subreso- 
nances are encountered as the system drifts can be deter- 
mined by considering the two resonant angles ip — w and 
4> — T^p. The time derivative '4> — w = jn~{j— l)np — w. 
As the planet migrates outward Up drops. For w > Wp 
the time derivative of the corotation resonant angle 
crosses zero first. For c > (corresponding to c < 
and a positive precession rate for the longitude of peri- 
apse), the corotation resonance is encountered first by a 
particle exterior to a planet as the planet migrates out- 
ward. This is what is expected for external resonances 
in a single planet system where the precession rate of the 
planet's longitude of periapse Wp = and dj > due to 
secular precession induced by the planet. 

We consider which resonance is encountered first for 
other drifting systems. Dust particles migrating inward 
and exterior to a planet would encounter the resonances 
in the same order, corotation resonance first as long as 
w > Wp. Here n is increasing, whereas for the planet 
migrating outward Up was decreasing. A particle located 
internal to a planet that is migrating inward would also 
encounter the resonances in the same order when m > 
Wp] (in this case one considers (j -\- l)np — jn — rij and 
(j -|- l)np — jn — Wp with Up increasing). 

We can compare the width of the corotation resonance 
to the separation between them. For the e-resonance, the 
libration width in A depends on the particle eccentricity 
or r. However since the corotation term does not depend 
on r, the corotation resonance width can more easily be 
estimated as AA ~ y/f. For 

|c| < Ve, 

the two resonances must overlap. Since c cx /i and 

y/i cx ep^^ii~^^^, strong resonances are likely to over- 
lap for small planet masses unless the planet eccentricity 
is extremely small. Because the two resonant terms dif- 
fer in sign, they have fixed points at different angles, and 
they are expected to interfere when overlapped even only 
slightly separated. 

We have found that a nonzero value of c does change 
the probability of capture. Figure |H1 shows numeri- 
cal measurements similar to that of Figure but for 
c = ±0.1. For c < (Figure IH^l) the corotation reso- 
nance is encountered after the e-resonance. The onset 
of the corotation resonance can knock the particle out 
of resonance, following capture into the e-resonance. For 
c > the corotation resonance is encountered first as 
the planet migrates (Figure Et*)- The capture proba- 
bilities are primarily modified at low drift rates where 
higher values of e are required to reduce the probability 
of capture. It is not obvious why this is the case. From 
individual integrations we note that the frequency c sets 
an oscillation period that is longer for smaller values of 
c. 'Holma n fc Murravl (|1996) found that c sets the Lya- 
punov time of the resonance. So for an overlapped system 
we might expect more highly chaotic behavior for larger 



values of c, particularly at low drift rates. Oddly higher 
values of c at low drift rates seem to stabilize the sys- 
tem, requiring higher values of the corotation resonance 
strength to kick the particle out of the e-resonance. 

For c < (as shown in Figure IHt) the corotation 
resonance is encountered after the e-resonance, conse- 
quently the corotation resonance is encountered after the 
e-resonance captures a particle. If the corotation reso- 
nance is strong it can knock the particle out of resonance. 
For larger separations, c ~ 1, temporary captures are fre- 
quent at large e. The extended low probability contours 
on the top end of Figure|Hb are in part due to temporary 
captures. 

4.1. Second order resonances 

When the resonance is second order it contains three 
subterms (see Equation IJ) . The first term (x F cos(^ — 
2w) and can be called (e.g., iMurrav & DermotHllQQfll) 

the e^-resonance since F cx at low eccentricity. The 
second term cx e^F^/^ cos^tp — w ~ Wp) and can be called 
an ee'-resonance. The third term cx cos(i/' — 2wp) can 

be called a corotation or e'^ resonance. Since the coro- 
tation term does not depend on F, its volume in phase 
space does not grow as the planet drifts and it should 
not be able to capture particles. However, as was true 
for the first order resonance, this resonant term can pre- 
vent the other resonant terms from capturing particles. 
When the corotation term is not large, both the e^- and 
ee'-resonances can capture particles. Since it is cx F^/^, 
the ee' term behaves like a first order resonant term, 
whereas the e^ term, cx F, is a second order term. In 
the previous section we found that first order resonances 
captured at a somewhat higher drift rate than the second 
order term and did not require as low initial momenta to 
exhibit a sharp transition between capture and no cap- 
ture. The critical drift rate for first order terms is cx //^/^ 
and for second order terms cx , a much steeper function 
of planet mass. Consequently it is possible that the ee' 
resonant term will capture particles and the e^-resonance 
will not capture particles even when the planet eccentric- 
ity is low. 

We first consider the situation where the e^-resonance 
is dominant. Taking k = 2 terms from Equation 10), we 
rescale the Hamiltonian as follows, 

K^K, V; r, 7) = A^ -h &5 A -f c^f + f cos(V' - 2w) (18) 
— C0S(V' — W — Wp) + cos(V' — 2Wp) 

where 



r = 


(52,0 

a 


-1 

F 




r = 


(52,0 |t 




c« = 




r^sij 


5n(a) 




r^sij 


5n(a) 




(524 ^2,0 






(^2, 211^2, oT 





The coefficients that can be adjusted are the strength 
of the ee'-resonance (set by f), the strength of the coro- 
tation resonance (set by e^), and the separation between 
the resonances (set by c^). Since the precession rate 



8 



c cx /i and the perturbation strength 62,0 fJ- we find 
that the resonance separation (c) does not depend on 
the planet mass. This impHes that the subterms could 
be well separated. This is different than the first order 



resonances that have 



C (X fj. 



1/3 



which implies that the 



resonance subterms are often overlapped. 

The ee'-resonance strength ^ oc fi^^^'^Cp. Since this 
depends on a negative power of /i, at low planet rnasses 
and at high planet eccentricities it is possible that ^ > 1. 
In this case the ee'-resonance could be dominant, and we 
would rescale momentum and time as we did for a first 
order resonance. In this case we could work with the 
Hamiltonian (equivalent to the previous one except for 
the rescaling) 

-A^ + b^A + c^f (19) 
+f cos{ip — w — zup) 
— XF cos(?/; — 2m) + cos('0 — 2mp) 



where 



>2.1 



-2/3 



r 

1/3^ 



X 



^'|<^2ar'/'|ar'/'sign(a) 



2/3|a|-i/3s 

2/3^1-1/3 



a\ "'"sign(a) 



= |'^2,o| 1(52,1 1 



(20) 

A comparison between Equation H19I) and Equation 
pnji shows that 

~e=x-' (21) 

We suspect_that the e^-resonance is more likely to cap- 
ture when ^ reaches a transitional value that we denote 
Straus- The ee'-resonance may be more important when 

X ^ ^trails- If C ^ ^trans then wc expcct sccoud order 
behavior {fip^crit oc /i^) otherwise we expect first order 
behavior {fip^crit oc /i^^'^). We note that the order the 
resonances are encountered is also important. Here we 
numerically measure a capture probability that does not 
specify which resonance captures. If the drift rate is suf- 
ficiently slow that the e^-resonance captures and this res- 
onance is reached first then this resonance will dominate 
the capture probability. 

As we did for the k = 1 resonances we have measured 
the capture probability for a range of coefRcients and 
drift rates. Figure |^ shows the capture probability for 
the Hamiltonian of Equation H18|) with varying drift rate, 
and ee'-resonance perturbation strength, ^, and with no 
corotation term; = 0. The resonances are not sepa- 
rated; = 0. On the lower left side of this plot, for 
weak ^ we see capture behavior consistent with the pure 
second order {k = 2) system discussed in the last sec- 
tion with only one resonant term. With the variables 
defined in this section (which differ by a factor of fc = 2 

from those defined in Equation [Tjl ^ =0.5. The 

crjt _ 

capture probability ceases to depend on ^ for ^ < 10~^. 
Consequently we can estimate a transition value 

-2 



valid for low initial particle momentum. For ^ > ^trans, 
the transition between capture and no capture occurs at 
faster drift rates and the transition drift rate is a function 
of ^. For ^ > ^trans the system behaves like a first order 
system, and since the perturbation strength depends on 
f we expect the transition drift rate to depend on f^/'^. 
This is consistent with the trend shown on the upper 
right in Figure |^. For large ^ and at high drift rates 
the system can fail to capture into the e^-resonance but 
can be captured into the ee'-resonance. In short, for 
£, ^ S,trans the systcm behaves like a second order system 
and tends to capture into the e^-resonance, however for 
C ^ £,tran8 the highly overlapped system behaves like 
a first order system and tends to capture into the ee'- 
resonance. 

FigureHjD shows the capture probability for the Hamil- 
tonian of Equation l(T!?)l with varying drift rate, varying 
e^-resonance perturbation strength, x with no coro- 
tation term; e-^ = 0. The resonances are not separated; 
Cx = 0. This figure extends Figure to ^ > 1 because 
X = At small x the system critical drift rate is 

independent of x ^'^'^ consistent with that measured in 



section 3 for a first order resonance with 



= 2. 



^tra 



10" 



(22) 



At this limit the system behaves like a first order reso- 
nance. The second order regime is not fully reached until 
£, ~ Sytrans which would correspoud to % > 20. However, 
for large values of ^ we see that the critical drift rate does 
begin to increases, consistent with a dependence of crit- 
ical drift rate on x (the second order term dominates). 
In Figure I^Id where x ~ 0.5 we see an extended region 
of drift rates corresponding to a regime of intermediate 
capture probability. We have inspected individual inte- 
grations from this region and seen temporary captures 
and widely varying or chaotic trajectories. The inter- 
mediate capture probability measured is probably due 
to interference between the two similarly sized and over- 
lapped e^- and ee'-resonances. 

Figures IHl shows integrations done for low initial parti- 
cle eccentricity (low f (tp)). We discuss what we expect 
would happen for higher initial particle momentum. The 
second order resonance is more strongly affected by the 
initial particle momentum. For larger r(to) the drop in 
probability at the bottom right on FigureEt would occur 
at a faster drift rate and the drop would be smoother. 
The contours would broaden and shift to the right, pri- 
marily on the bottom of Figure The resonance be- 
haves like a second order resonance at the top right of 
Figure So we expect a similar broadening and shift- 
ing to the right of the contours at the top of Figure EJa. 

We now investigate the role of the corotation or e'^ 
resonance term. When the e^ term is not important (as 
we have found ^ > itrans), the system is identical to 
that studied in the previous section for fc = 1 with two 
terms fEauation ll6|l . However we can also study the sys- 
tem for Equation 1181 as a function of and with f = 
so the corotation term is strong but the ee'-resonance is 
weak. Figure 1101 shows the capture probability for the 
Hamiltonian of ^| with varying drift rate, and 
with £^ — — 0. We find here that even moderate values 
of the corotation resonance strength e '--^ 0.1 can signif- 
icantly reduce the capture probability. The corotation 
resonance can reduce the probability of capture for both 
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first and second order resonances. 

4.1.1. Separated second order resonances 

As was true for the first order resonances we expect the 
capture probabiHty to depend on the resonant term sep- 
arations and order that the resonances are encountered. 
Figure^Jshows the effect of changing the resonance sep- 
aration when the e^-resonance dominates. Figure [TTk 
shows the case with widely separated resonances when 
the e^-resonance is encountered first. We see that the 
transition value of S,trans is higher than when there is no 
separation {c — 0). When the e^-resonance is encoun- 
tered afterward the transition value of ^trans is lower. 
The ee'-resonance interferes with the capture into the e^- 
resonance to a higher degree when this resonance is en- 
countered earlier. Figure ^2 shows the effect of changing 
the resonance separation when the ee'-resonance dom- 
inates. We find that the region of intermediate capture 
probability at x ~ 0.5 is smaller when the resonances are 
separated than when c = (Figure Et). 

5. APPLICATIONS 

In this section we apply what we have learned above 
to two systems involving capture into the 2:1 resonance. 
When the 2:1 resonance is exterior (and capture particles 
as a planet migrates outward) the e-resonance strength 
is reduced because of the indirect term. This reduces the 
critical migration rate compared to that for other first or- 
der resonances. Because the corotation resonance is not 
affected by the indirect term, it is comparably strong. 
Consequently even a small planet eccentricity can reduce 
the capture probability. When the 2:1 resonance is in- 
terior (and can capture for a planet migrating inward) 
the indirect term reduces the strength of the corotation 
resonance instead of the e-resonance. In this case the 2:1 
resonance is strong and can capture at fairly high mi- 
gration rates. We consider two situations, the capture 
of twotinos into the 2:1 resonance by Neptune migrating 
outward, and the capture of an inner extra solar planet 
into the 2:1 or 3:1 resonance by an inward migrating 
planet exterior to it. 

5.1. The capture of twotinos in the Kuiper Belt via 
Neptune 's migration 

In this section we consider the capture of Kuiper belt 
objects into the 2:1 resonance by an outward migrating 
Neptune. We see from Table ^ that the 2:1 external 
resonance has exceptionally large values of e compared 
to the 3:2 and 4:3 resonances. This is because Si^ is 
small due to the addition of the indirect term. 

We first consider the critical migration rate allowing 
capture and compare the one we predict here with that 
found from numerical studies. Migration rates are often 
given in terms of the time it takes to cross the range 
of radius covered during the entire migration . This is 
typically a few AU for Neptu ne's migration l|Ida et al.l 
l200nt IChiang fc ,ToTdai3 l200l . To compare migration 
rates to the critical one estimated above, we must first 
convert rates into our system of units. For GMq — 1 
and radii in units of Neptune's semi-major axis, ajv, we 
multiply timescales by yja^j /GMq — 26.1 years. A mi- 
gration rate of a few AU in 10^ years corresponds to 
hp ^ 2.6 X 10^'' f ^ IP V Since oc a^'^^^ this corre- 



sponds to Tip ~ 3.9 X 10 y ^ ^" ^'^ j . The critical planet 
migration rate (listed in Table ^ for Neptune's 2:1 reso- 
nance IS Up^crit ~ 

0.54//3 ^ i Q X 10-6. We find that 
J^^O.J . (23) 

'^p.crit \ ^migrate J 

IChiang fc Jorda n (2002) found that 2:1 resonance cap- 
tured at the 50% level for t,nigrate = 10^ years but was 
much less efficient, capturing only 15% of particles at 
faster migration rates of tmigrate ^ 10^ years. The same 
rise in c apture probabili ty at tmigrate ~ 10^ years was 
seen by llda et all l)200Cl() . The sharp drop in capture 
probability is consistent with our predicted limit for the 
critical migration rate. We find that we ca n account 
for th e trends seen in the nume rical studies of llda et alJ 
l)200nt) : IChiang fc ,TordanH20r)^ and confirm the theoret- 
ical explanation of IfMedlandf |200lD . We note that the 
transition from a probability of 50% to 15% occurs over a 
fairly large range of drift rates. In our toy model we could 
account for such a smooth transition with initial parti- 
cle eccentricity near eum- However because the capture 
probability drops steeply for eo > eum this explanation 
would require fine tuning of the initial particle eccentric- 
ity distribution. We note that it is impossible to zero the 
eccentricity of a particle in a simulation because of other 
perturbations. Also, because we have dropped most co- 
sine terms in Equation Q , We have neglected these other 
perturbations in our Hamiltonian model. Consequently 
it is difficult for us to compare the initial eccentricity 
distribution of a simulation to the distribution in our 
momentum F. 

We now consider the role of the corotation resonance. 
From table ^ we find e = 6.6fi~^^^ep. For Neptune 
^ = 5.1 X 10~^ and we find that e « ISOep. For Nep- 
tune's current eccentricity ep ~ 0.008 this places e ^ 1.4. 
This is somewhat above the critical corotation strength 
value allowing capture into the 2:1 resonance according 
to figure[7|when the resonances are on top of each other. 
We need to consider the separation between the reso- 
nances; c = 1.18/x^/'^ — 0.04. However this is the sep- 
aration only if Neptune's longitude of perihelion does 
not precess. Neptune's precession frequency is largely 
due to the solar system's eighth eigenvector that dom- 
inates this planet's se cular motions (Nobili et al. 198(| 
lApplegate et al.lll986(l . Neptune's precession rate due to 
other planets is a few times larger than that it induces 
on objects in its 2:1 resonance. Consequently c could be 
larger a factor of a few and either positive or negative de- 
pending upon the secular motion of the planet when the 
migration took place. Smaller values of planet eccentric- 
ity would allow the 2:1 to capture whereas larger values 
would tend to make it more difficult. Negative values of 
c would allow the resonance to be in the temporary cap- 
ture regime shown in Figure|SjL whereas positive values of 
c would allow capture at Neptune's current eccentricity, 
provided the migration was slow (Figure IHb). It is in- 
teresting to find that Neptune's eccentricity is very near 
the critical value that would make this resonance fail to 
capture. This large value of e could contribute to the 
intermediate capture probability seen in simulations and 
the moderate range of drift rates where this intermediate 
capture occurs. 
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5.2. Capture into the 2:1 and 3:1 Resonances of 
Multiple Extrasolar Planet systems 

Three extrasolar multi-planet systems have two plan- 
ets locked i n the 2:1 resonance , GL876, HD 82943, and 
HD128311 ijMarcv et al.ll2005j) . In each case the outer 
planet is more massive than the inner one. The masses 
of the outer planet are 1.9, 1.6 and 3.2 2 Mj (Jupiter 
masses), respectively. We assume that an outer planet 
has migrated inward and captured the interior and lower 
mass plane t into the 2:1 resonance (e.g., as explored pre- 
viously hv IKlev et a.LlIM 1200^ iMoorhead fc Ad'a^^^^ 
120051) . The coefficients for this situation are listed in the 
appendix and in Tabled For an internal 2:1 mean mo- 
tion resonance fip^crit = 22.7/x'*/'^. We relate the critical 
mean motion drift rate to a critical semi-major axis drift 
rate (with a factor of 2/3) and restore the units. To cap- 
ture an internal planet into the 2:1 resonance a planet 
must have a migration rate slower than 

«P < 15a^^/^ ( ^] ■ (24) 



Using the period of the planet's orbit, P — 2t:\J -^j^, 

we can relate the critical migration rate to a timescale, 
T — Op/ hp, finding 

T2:l > 0.4^-4/3p (25) 

for the 2:1 resonance For a 2Mj mass planet, we find a 
migration timescale of longer than 1600 orbital periods 
is required for the 2:1 resonance to capture. This limit 
is consistent with the timescales a dopte d for migration 
in the simulations by IKlev et all 1)20041) . We can con- 
sider the eccentricity limit (taking the value from Table 
Q]l eura ~ 1.5/i^/^ = 0.2. This implies that the initial ec- 
centricity of the inner planet (as long as it was below 0.2) 
would probably not limit the capture probability into the 
2:1 resonance. 

We now consider capture into the 3:1 resonance. The 
55 CnC system has two planets locked in the 3:1 res- 
onance with outer planet 0.2Mj and inner planet with 
0.8Mj. Even though the inner planet is more massive 
we consider capture into an internal resonance because 
the 3:1 internal resonance is stronger than the external 
one. This is a result of the contributions from the indi- 
rect terms. Using values given in the appendix in Table 
|21for an internal 3:1 resonance we find hp^crit = 0.6/i^ for 

the resonance and fip crit = 4:lfi'^/'^ef/^ for the ee' res- 
onance, in both cases for low initial particle eccentricity. 
This corresponds to a migration timescale 

r3:i,e2 > 15//-2p (26) 

for the resonance and 

T3:l,ee' > 0.2^,-^/%^/' P (27) 

for the ee' resonance. However the migration rate is 
less restrictive for the e^ resonance if the inner planet 
has a moderate eccentricity. The limiting eccentricity is 
eu^ - 0.2^1/2. For fi = 0.002, we find e^™ = 0.003. 
Consequently the inner planet is likely to have e > enm- 
In this case the migration timescale must be modified 
by the factor given in Equation p5(l and the limiting 
migration timescale would be ^ 10 times smaller or 

T3:l,ee' > 0.02M-4/3g^4/3p_ (38) 



For ^ — 0.0002 to capture we find the migration 
timescale must be longer than a few times 10^ orbits 

for the e^ resonance and 20000ep ^^'^ orbits for the ee' 
resonance. For a moderate planet eccentricity of 0.2 this 
corresponds to a timescale of 2 x 10^ orbits. 

We have found here that capture into the 2:1 resonance 
by the multiple planet extra solar systems does not re- 
quire a slow migration rate but capture into the 3:1 res- 
onance does. For the 3:1 resonance faster migrations are 
allowed for the ee' resonance than the e^ resonance. The 
limiting timescale is decreased if the planet eccentricities 
are not low prior to resonance capture. 

6. SUMMARY AND DISCUSSION 

In this paper we have explored the problem of res- 
onance capture for mean motion resonances at fast or 
non-adiabatic drift rates. We first studied the first and 
second order time dependent Hamiltonian system with 
one resonant term. We find that for sufficiently low ini- 
tial particle momentum (or eccentricity), the transition 
between resonance capture and no capture is sharp, oc- 
curring over a narrow range in drift rate. We give an ex- 
pression fEauation ll4(l which makes it possible to predict 
the critical planetary migration rate (above which there 
is no capture) for first and second order mean motion res- 
onances in the general restricted three body problem in 
the limit of low initial particle eccentricity. Expressions 
are given in the appendix for coefficients which allow one 
to estimate the critical drift rate for any first or second 
order mean motion resonance. Coefficients are evaluated 
for strong resonances and listed in Tables [H and El This 
genera lizes upon previous analytical work by Fricdlan^ 
()2001l) and provides a theoretical explanation for critical 
drift rates measured numerically and their depe ndence 
on planet mass fe.g.. llda et al.ll200nt rW^ '2()0?). 

We have numerically measured the probability of cap- 
ture as a function of initial particle eccentricity. We find 
that the transition between resonance capture and no 
capture is smoother, occurring over a larger range in drift 
rate, for initial particle eccentricity of order the limit- 
ing value ensuring capture in the adiabatic limit, eum- 
The drift rate at which the capture probability is half is 
not strongly dependent on the initial particle eccentricity 
for first order resonances, and the probability of capture 
drops rapidly for initial particle eccentricities exceeding 
the limiting value, eum- For second order resonances, we 
find that the drift rate at which the capture probabil- 
ity is half is higher when the initial particle eccentricity 
is higher. Equation 1)15(1 can be used to estimate the 
half probability drift rate for initial eccentricities below 
10e;im. At eo ^ SOeum the capture probability drops 
below 1/2 at all drift rates. 

In the limit of low initial particle eccentricity, we have 
considered the case of resonances containing multiple 
subterms. A first order resonance fails to capture when 
the corotation resonance has unitless strength e>l. As 
this coefficient depends on planet eccentricity, migrat- 
ing, eccentric, low mass planets could have first order 
resonances that fail to capture particles for this reason. 
A regime of intermediate capture probability also ex- 
ists at high drift rates and large corotation perturbation 
strength. We have found that the resonance separation, 
and order of encounter, affects the capture probability, 
primarily when the corotation resonance is strong. This 
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implies that the capture probabihty is dependent upon 
the precession rates of the longitude of periapse of both 
particle and planet. 

Second order resonances contain three subterms. As 
was true for the first order resonances, the corotation 
resonance fails to capture particles but can prevent the 
other resonant terms from capturing particles if scale free 
parameters £5 ^ 1 or > 1. This implies that above a 
certain planet eccentricity, second order resonances fail 
to capture particles. Below this planet eccentricity the 
e^- and ee'-resonances can capture particles. When our 
coefficient ^ < 10~^ the e^-resonance will capture parti- 
cles (providing the drift rate is sufficiently slow) and the 
capture behavior is second order. For x < 0.1 the ee'- 
resonance will capture particles and the behavior is first 
order. For nonzero planet eccentricity, the ee' subreso- 
nance (which behaves like a first order resonance) may 
more easily capture particles at faster drift rates than 
the resonance. For x or ^ of order 1, a regime of in- 
termediate capture probability exists at high drift rates. 
For second order resonances the subresonance separation, 
and order of encounter also affects the capture probabil- 
ity. 

A number of effects have been proposed to account 
for reduction in capture probabilities compared to those 
predicted via adiabatic theory, e.g., Zhou ct al. (2002) 
showed that stochastic or jumpy migration would allow 
particles to escape resonances. Here we have shown that 
rich dynamics in the non-adiabatic limit allows particles 
to escape resonance capture. We have shown that corota- 
tion terms can reduce the capture probability. For second 
order resonances, resonant subterms can interfere, again 
producing a regime of intermediate capture probability. 
For first order resonances, the half probability drift rates 
are not strongly dependent on the initial particle eccen- 
tricity, and the probability of capture drops rapidly above 
a limiting initial eccentricity. However for second order 
resonances the half probability drift rate is higher for 
initial particle eccentricity near the limiting value. Con- 
sequently we expect that second order resonances should 
have larger regimes of intermediate capture probability 
in range of drift rate and initial particle eccentricity. 

We have applied our understanding to the problem of 
capturing twotinos via Neptune's migration. We find 
that the eccentricity of Neptune is sufficiently high that 
the 2:1 resonance could fail to capture particles. Cer- 



tainly if Neptune's eccentricity were any higher during 
migration its 2:1 resonance would not have captured par- 
ticles efficiently. It is interesting to find that Neptune's 
eccentricity is very near the critical value that would 
make this resonance fail to capture particles. 

We have applied our framework toward predicting min- 
inum migration timescales allowing extra solar multiple 
planet systems to capture into the 2:1 or 3:1 resonances. 
We find that a migration timescale of greater than a few 
thousand orbital periods is required to allow capture into 
the 2:1 resonance for three systems. However a much 
longer timescale, ^ 10^, orbital periods is required to 
allow capture into the 3:1 resonance for the 55 Cnc plan- 
etary system. The migration timescale can be somewhat 
reduced if the planets are on moderately eccentric orbits 
subsequent to migration. 

In this work we have extended the theory of resonant 
capture for drifting Hamiltonian systems to the non- 
adiabatic limit and to systems with multiple resonant 
subterms. We have provided a theoretical framework 
to predict resonance capture probabilities. However this 
framework is based on numerical integration of a sim- 
plistic two-dimensional Hamiltonian model and so may 
not accurately represent the full complex dynamical sys- 
tems. Direct numerical integration of these systems must 
be carried out to test the validity and accuracy of the 
expressions given in this paper. The exploration done 
here could also in future be extended via numeric al study 
of mo dified quasi Hamiltonian toy models (e.s.. \GomeA 
((ini3)) to better cover systems with drift induced by 
non-conservative forces such as gas drag or Foynting- 
Robertson drag. This work could also be extended to 
cover motions out of the plane and high eccentricity sys- 
tems. 
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APPENDIX 

COEFFICIENTS FOR INTERNAL AND EXTERNAL RESONANCES 

External resonances 

For external first order resonances, 

K{A,tP;T,j)^aA^ + bA + cT (Al) 

+'^i,or^^^ cos — w) + di^i cos (V' — Wp) 

with 

Si,i^~Hepaf27. (A2) 
Coefficients a,b, and c are given in Equations Q and Q. The fi are functions of the Laplace coefficients and are 
evaluated at a with index j using expressions from the appendix bv lMurrav fc Dermo tt (1999). 
The above expressions only include direct terms. For the 2:1 resonance the indirect term contributes and 



5i.o(2 : 1) = -fiV2a^/^ L/31 - ^ 



(A3) 
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The near cancellation of the direct and indirect terms makes second order terms important for the 2:1 resonance 
l)Friedlandll2n01tlMurrav-Clav fc Chiangi r2005'l. It may be useful to recall the maximum initial particle momentum or 
initial eccentricity, eum ensuring capture in the adiabatic limit. This corresponds to 



where Tq,;™ = 3/2. 
For external second order resonances or fc = 2, 



1/4 



<5l,0 



1/3 



where 



K{K, ij; r, 7) = aA^ + 6A + cF 

+^2,oF cos (Vj — 2n7) 

+52,lF^/^ cos {-tjj — W — TUp) 

+62,2 cos (-0 — 2wp) 



h,2^ -f^elaf 45- 



(A4) 



(A5) 



(A6) 



Coefficients a,b and c are given in Equations 10) and For the 3:1 resonance the indirect term contributes and 

3 



^2,0 



8a 



1/2 



^2.0 (3 : 1) = -^2^1/2 l^^f^ 
For second order resonances capturing into the subterm 

e/im,{ = ^2fo,;i,„a^^'* 

where the critical scale free momentum Fo^hm = 1/8. For those capturing into the ee' subterm 



where the critical scale free momentum Fo^iim = 3/2 



2Fo^;imQ; 



1/4 



02,1 



1/3 



(A7) 



(A8) 



(A9) 



Internal resonances 

To make our theory appropriate for internal resonances (external perturber) we would consider j : j + k resonances 
and change the coefhcients to 







a = 




b = 


-{j + k){np-l) 


c = 


-m2/2«-^/' 


<5i,o = 


-MV2a-i/V27 


^1,1 = 




62fl = 




^2,1 = 


-[iep\Pia~^l^ Jasi 


<52,2 = 





(AlO) 

where a = ajap and we have used the approximation ~ 2T / L ^ 2Ta^^^^. The c t erm describes secular prece ssion 
of the longitude of periapse and depends on the function /2 given in the appendix bv lMurrav fc DermottI |1999) and 
is evaluated at a with index j = 0. For internal resonances the oth er fi functions are evaluated at a with index j + k 
using expressions from the appendix bv lMurrav fc DermottI ljl999t) . For the 2:1 and 3:1 resonances the indirect terms 
contribute and 

<5i,i(2 : 1) = -^ep(/3i-2a) (All) 



^2,2(3 : 1) = -fJ-el 



f 27 

753 - 



(A12) 



The maximum initial particle eccentricity ensuring capture in the adiabatic limit for first order resonances 

oi,o 



-1/4 



where TQ^um = 3/2. For second order resonances capturing into the subterm 



2fo. 



-1/4 



S2A 



1/2 



where the critical scale free momentum Tq^, 



1/8. For those capturing into the ee' subterm 



eii 



2ro, 



-1/4 



1/3 
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(A13) 



(A14) 



(A15) 



where the critical scale free momentum Fo^jim = 3/2. The drift rate for a capture probability of 1/2 given in Equation 
IjlSI) must be modified for second order internal resonances; 



^p,l/2| 



0.5(j - 2)5' 1 



02,0 



0.25 



(A16) 



This expression is valid for initial particle eccentricity eo ^ lOeiim,^- 
The coefficients for strong internal and external resonances are listed in Tables ^ and |21 

MIGRATION OF DUST VIA POYNTING-ROBERTSON DRAG 

In this paper we have considered an varying Hamiltonian system. However there may be some similarities between 
this system and the slowly drifting dissipative systems. We add relations that allow the reader to approximately predict 
the critical drift rates for dust spiraling inward under Poynting-Robertson drag. In the case of Poynting Robertson 
drag dust particles in a circular orbit decay on a timescale proportional (3~^ where (3 is the ratio of radiation to 
gravitational (from the star) forces. It is convenient to write 



0.2 



Ah 



(Bl) 



where is the radius of the particle in iim and is the luminosity of the star ijSicardv et alJlT993j) . The drag force 



leads to a slow increase in the mean motion 



3ai/2/3 



ci 



(B2) 



where c; is the speed of light in units of the planet's velocity or divided by ^ GlvUjap. The value of our coefficient b is 

not important, as long as it passes through zero on resonance. However, its drift rate or h is important. At resonance 
jn = (j — k)np and we can relate the drift rate of the particle spiraling inward to a system of a planet migrating 
outward considered in the previous sections. We replace hp with n, finding an effective drift coefficient 



The rescaled speed of light 



Consequently we can write 



ci 



10* 



b = 



ci 

-1/2 



\1AUJ 



1/2 



6 = 0.6 X 10" 



1/2. 



Me 



-1/2 



-1/2 



(B3) 



(B4) 



(B5) 



The above relation can be used to approximately determine the minimum size particles that can be captured into 
resonances using the formulation presented in previous sections of this paper. 
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Fig. 2. — Capture probabilities for a Hamiltonian system Equation JSJ with one first order (fc = 1) resonant term as a function of drift 



rate 



db 
dT 



and initial momentum r(io). Note that low r(io) corresponds to low initial eccentricity. The thick and thin solid lines shows 

the capture probability for f(io) = 10~* and lO"'^ respectively. The squares and stars show the capture probability for f(to) = 1 and 2.3 
respectively. In the adiabatic limit for r(to) < 1.5 the capture probability is 1. For r(to) = 2.3 the capture probability is intermediate for 
low drift rates and approaches a constant value as the system becomes more adiabatic. For initial momentum low (10~*) the transition 
between 100% capture and 0% capture is extremely sharp. We find that if the initial momentum is ~ 1 then there is a regime or a range 
of drift rates where the capture probability is intermediate. With a change of scale, all first order resonances can be put in the form of 
Equation jSJ. Consequently the probabilities shown here can be used to estimate the capture probability in the non-adiabatic limit for any 
migrating first order resonance. 
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log(d b/d t) 

Fig. 3. — Similar to Figure ISl except for the second order resonances (Equation |S] with k = 2). The thick, intermediate and thin solid 
lines show the capture probability for r(to) = 10~®,10~* and 10~^ respectively. The squares and stars show the capture probability for 
r(to) = 0-1 S'lid 1.0 respectively. Capture is ensured in the adiabatic limit for r(to) < 1/8. As was true for the first order resonances the 
transition between capture and no capture is steeper (covering a narrower range in drift rate) for low initial momentum. For second order 
resonances the drift rate for a capture probability of 0.5 depends on the initial momentum. 
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Fig. 4. — a) The drift rate at which the probability of capture is 1/2 (stars) and 1/4 (squares) as a function of initial momentum (r(io) 
) for first order resonances {k = 1). The x-axis shows log]^Qr(to). The j/-axis shows logj^g |^|- b) Width of the drop in probability as a 
function of initial momentum. 
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a) The drift rate at which the probabiUty of capture is 1/2 (stars) and 1/4 (squares) as a function of initial momentum for 
first second order resonances (fc = 2). The x-axis shows IogiQr(4o). The j/-axis shows log^Q |^|- The dotted line shows the function given 
in Equation 1121 . b) Width of the drop in probability as a function of initial momentum. 
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Fig. 7. — The capture probability for a first order resonance (fc = 1) is shown as a function of drift rate |^| and strength of secondary 
resonance e. Contours are shown at probabiUties of 0.15, 0.35, 0.55, 0.75 and 0.95. The numbers 0.0 and 1.0 on the top right and lower left, 
respectively are placed to make it clear where the probability of capture is near zero and near 1. Equation 1161 was numerically integrated 
with resonance separation 5 = 0, and low initial momentum T{to) = 10~*. Each probability was measured from 100 different trials 
with randomly chosen initial angles. The x-axis shows logjg |^| and the y-axis shows logj^Q(e). For low values of corotation perturbation 
strength, e, the transition between capture and no capture happens at the critical value of estimated for the case of a single perturbation. 
Near the critical drift rate the additional resonant perturbation can cause a moderately large region with an intermediate probability of 
capture. At lower drift rates, the corotation resonance prevents capture for e ^ 1. To the upper right on this plot (large drift rates and 
large e values), temporary capture can take place, making it more difficult to measure absolute capture probabilities. 
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Fig. 8. — a) Same as Fig |7| except the resonance separation c = —0.1. The frequency of the corotation term is shifted so that this 
resonance is encountered after the e-resonance. The onset of the corotation resonance can kick the particle out of resonance, b) Same as 
Fig|7|except the resonance separation c = 0.1. The frequency of the corotation term is shifted so that this resonance is encountered first. 
At lower drift rates larger e is required to prevent captures. These figures show that the subresonance separation can influence the capture 
probability. The subresonance separation is set by the difference between the planet and particle's precession rate of longitude of periapse. 
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log(db/dT) log(db/dT) 

Fig. 9. — a) The capture probability for a second order resonance (k = 2) as a function of drift rate and ee'-resonance strength. This 
figure is similar to Fig. |7| The Equation 1181 was numerically integrated with resonance separation = 0, no corotation term, = 0, and 

initial momentum f (4o) = 10^^. Drift rate and f were varied. The x-axis is logjQ |-^| and the i/-axis is logjpd^). For low ee'-resonance 
strength or f < 10""^, the capture probability drops at a drift rate consistent with the critical value measured in section 3 for a single 
second order resonance fEguation llli . For larger values of ee' resonant term, ^, the critical drift rate increases, depending upon ^. In this 
case the system fails to captu re in to the second order e^-resonance and instead captures into ee'-resonance which behaves like a first order 
resonance, b) The Equation 1191 was numerically integrated with = 0, e;^ = 0, and initial momentum r(to) = 10~^. Drift rate and x 

were varied. The x-axis is logj^Q ^'^'1 ^^'^ j/-axis is logiQ(x). Because x = g"^/^ this figure covers large values of extending past the 

top of a). For low e'^-resonancc strength or § < 10 ~^, the capture probability drops at a drift rate consistent with that predicted in section 
3 for a single first order resonance; ('Eguation lili . For x ~ 0-5 there is a regime of drift rates with intermediate capture probability. 




-1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 

log(db/dx) 



Fig. 10. — The capture probability for a second order resonance (fc = 2) as a function of drift rate and corotation resonance strength. 
This figure is similar to Fig. |7| The Equation 1181 was numerically integrated with = 0, § = 0, and initial momentum T{to) = 10~®. 

Drift rate and ej were varied. The x-axis is log^Q |-^| and the y-axis is log^Q(e^). At low values of the capture probability drops at 
a value consistent with that predicted in section 3 for a single second order resonance; Equation 1111 . At ~ 1 the corotation terms 
prevents capture into the e^-resonance. 
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Fig. 11. — Separated second order resonances, a) This Figure is similar to Figure|5Ji except 5^ = —1. The ee'-resonance is encountered 
after the e^-resonance. b) This Figure is similar to Figure [HJ. except = 1. The ee'-resonance is encountered before the e^-resonance. 




Fig. 12. — Separated second order resonances, a) This Figure is similar to Figure|5J) except c-^ = —0.5. The ee'-resonance is encountered 
before the e-^-resonance. b) This Figure is similar to Figure|2l3 except = 0.5. The ee'-resonance is encountered after the e^-resonance. 
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TABLE 1 

Coefficients for k = 1 resonances 





Exterior 


resonances 


Interior 


resonances 


j : k 


4:3 


3:2 


2:1 


4:3 


3:2 


2:1 


a 


0.825 


0.763 


0.630 


0.825 


0.763 


0.630 


a 


-16.35 


-7.86 


-2.38 


-19.81 - 


-10.30 


-3.78 


c/m 


-4.14 


-2.01 


-0.61 


-5.02 


-2.64 


-0.98 


<5i,o/ai 


-3.41 


-2.31 


-0.24 


-4.21 


-3.06 


-1.89 


Si,i/ilJ.ep) 


2.35 


1.55 


0.75 


3.07 


2.29 


0.31 


"'P,crzt/ H' 


198.5 


48.4 


0.54 


398.7 


126.4 


22.7 


e7(^-l/3ep) 


1.16 


1.00 


6.59 


1.22 


1.21 


0.21 


c//il/3 


0.72 


0.58 


1.18 


0.71 


0.58 


0.41 




0.98 


1.08 


0.72 


1.08 


1.24 


1.54 



Note. — Equations and jlj were used to c alcul ate th e coefficients 
a and c (for external resonances). Equations iA2l and IA3I were used to 
estimate the perturbation strengths for the 4:3 and 3:2 resonances and 
the 2:1 resonance, respectively, for external resonances. Coefficients for 
the internal resonances are g iven in the appendix. The critical drift rate 
is calculated fro m E quation 1141 . The coefficients e and 5 are calculated 
from Equation 1171 . We have calculated c using only the secular term 
from one planet and assumed that rijp = 0. The critical eccent ricity 
ensuring capture in the adiabatic limit is calculated using Equation 1A41 . 



TABLE 2 
Coefficients for k = 2 resonances 



j : k 


Exterior 


resonances 


Interior 


resonances 


3:1 


5:3 


3:1 


5:3 


a 


0.481 


0.711 


0.481 


0.711 


a 


-3.12 


-18.98 


-6.49 


-26.68 


c/n 


-0.20 


-1.22 


-0.41 


-1.72 


'52,o//i 


0.24 


6.82 


0.63 


3.47 


<52,l/(Aiep) 


-1.25 


-8.00 


-3.76 


-13.33 




0.10 


1.04 


0.36 


5.69 




18.56 


1.95 


19.11 


10.61 


6-5/(M-ie2) 


5.59 


0.43 


5.93 


12.56 




0.81 


0.18 


0.65 


0.50 




0.029 


69.86 


0.60 


30.19 




0.40 


0.95 


0.19 


0.20 




0.114 


0.174 


0.116 


0.538 


C-x/{/.V3,;2/3)^^^ 


0.12 


0.11 


0.09 


0.10 




5.77 


683 


40.66 


1694 




1.06 


1.19 


1.73 


1.50 



Note. — Equations j2J and jlj were used to calculate a and c. The 
critical drift rates hp crit,^ ^'nd np^crit,x ^-re those calculated for low ^ 
and X respectively, »^p_cr^t,? given in the limit of low initial particle 
eccentricity. For particle initial eccentricity near the limit ensuring cap- 
ture in the adiabatic limit, the values given here for np ,.^it^^ should be 
multiplied by 10. The coefficients e^, and are calculated from Equa- 
tion 1191 . The coefficients i ^i^d ^^e calculated from Equation 1201 . 
Expressions from the appendix were used to calculate the <5 coefficients 
and critical eccentricities. We have calculated c using only the secular 
term from one planet and assumed that vop = 0. 



